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This article presents the exact potential corresponding to confining fields inside 
a linear rf quadrupole particle trap of finite length. The analytic expression for 
the trapping potential is derived by introducing a linear trap employing a relatively 
simple cylindrical geometry and solving Laplace’s equation for the trap electrodes. 
The finite length of linear traps results in field distortion near the trap ends. An 
exact analytic determination of the fields is useful because the profile of the trapped 
ion cloud is highly dependent on the fields confining it. It is shown that near the 
ends of the trap , the effective potential arising from the rf fields acts to propel 
particles out of the trap , and further, that the addition of a dc bias generates an 
inhomogeneous field in the trap that influences the particles both perpendicularly 
to and along the trap’s long axis. 

I. Introduction 

Trapped 199 IIg + ion standards are presently the most 
stable frequency standard developed for averaging times 
> 10 4 seconds. Trapped ions can be stored for long pe- 
riods, essentially isolated from the outside environment. 

This arrangement serves to reduce the perturbation to the 
ions’ atomic-energy levels and hence minimize frequency 
fluctuation. 

Currently, the Time and Frequency Systems Research 
Group is developing a prototype linear ion trap consisting 
of four parallel cylindrical rod electrodes enclosed between 
two end electrodes (Fig. 1). This trap is preferable over the 
conventional spherical trap in frequency-standard applica- 
tions because of its 20-times-larger ion storage capacity. 

The number of ions determines an important trade-off pa- 


rameter between the signal- to-noise ratio and the second- 
order Doppler shift arising from the motion of the ions 
caused by the trapping field. The geometry of the linear 
trap improves clock peformance by allowing an increase 
of the number of trapped ions without degradation of the 
frequency stability [1]. 

Particle motion is directly affected by the electric fields 
generated within the trap. A time-varying quadrupole 
electric field is produced in the linear rf trap by an ac volt- 
age applied to the four rods such that any rod’s nearest 
neighbor has the opposite polarity. The rf field provides 
the trapping force necessary to confine particles in the ra- 
dial direction (see Fig. 1). In order to prevent trapped par- 
ticles from escaping out the ends of the trap, a dc voltage is 
applied to the two “endcaps ” producing an exponentially 
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decaying field directed inward along the the axis of cylin- 
drical symmetry of the trap (i.e, along the z axis). The 
linear trap is also capable of containing larger macropar- 
ticles (such as alumina or other substances) [2]. In this 
case, however, static electric fields must be incorporated 
in addition to the rf and endcap fields to offset gravita- 
tional forces on the macroparticles, which are comparable 
to the trapping forces. If the trap is oriented such that 
the trap axis is horizontal with the rods at the corners of 
a square when viewed down the trap axis, a dc voltage 
can be applied to the top two rods of the quadrupole to 
provide this offset. 

RF fields in the radial direction of a linear trap have 
been analyzed extensively by assuming the trap to be in- 
finitely long [3]. In actuality, however, the traps are of 
finite length, and thus both static and time- varying fields 
emanating from the rod electrodes will be directed lon- 
gitudinally along the axis of the trap, as well as radially. 
Near the ends of the trap, field distortion occurs. In an ion 
trap, this can subsequently aflect the shape of the trapped 
ion cloud and thus ion motion. To optimize clock per- 
formance, an exact determination of the fields generated 
within the trap, including the “end effects,” is necessary. 
This then allows for the calculation of the surface shape of 
the ion plasma contained within the trap. In this article, 
an analytic model is developed for one implementation of 
a linear trap to determine static fields produced by the 
end electrodes and time- varying electric fields everywhere 
inside a finite-length linear trap with cylindrical geome- 
try. The model consists of a hollow cylinder, partitioned 
lengthwise into four sections (analogous to the four rods) 
and two endcaps. The simple geometry of the cylinder 
allows the calculation of the fields within the cylinder ac- 
curately (Fig. 2). 

Although consideration of a dc offset to counterbalance 
gravitational forces is unnecessary when modeling a linear 
trap designed for the storage of ions, it cannot be neglected 
when discussing a trap that contains macroparticles. The 
primary purpose of applying the dc offset is to counter- 
balance gravitational forces on macroparticles. However, 
if the bias is added to the rod electrodes, electric fields 
will be additionally generated along the axis of the trap in 
a manner that acts to longitudinally localize the particles 
in the trap center or to propel them out the ends of the 
trap, depending on the polarity of the dc offset. A second 
model is presented in this article to determine the effects of 
the fields generated by a dc voltage applied to two rods in 
the linear trap. Because the focus of this model is solely 
to demonstrate the effects of applying a dc offset to the 
trap, the contribution from the rf trapping fields and the 
fields generated by the end electrodes are neglected, al- 


though the complete determination of the potential inside 
a cylindrical macroparticle trap is possible by combining 
the potentials derived from the two models. 

First presented is the conventional trapping theory 
for rf quadrupoles previously developed in the literature. 
Then, the potential everywhere inside a linear ion trap 
with cylindrical geometry is derived by modeling the trap 
after a four-sectioned hollow cylinder of finite length en- 
closed by two endcaps. In Section IV, the potential is 
determined inside a second model that consists of a hol- 
low cylinder split along its length into two halves. The 
top half is at a voltage V 0 and the bottom half and end- 
caps are grounded. This model accounts for the effect of 
a dc-offset potential applied to the rod electrodes in a lin- 
ear macroparticle trap. The conclusion summarizes the 
results of the two models. 


II. Trapping Theory 

The linear trap typically employs two dc-biased “end- 
cap” electrodes to confine the particles longitudinally, and 
four rods that produce a time-varying electric field to con- 
tain them in the transverse or radial direction. The volt- 
age applied to the rods is such that at a given instant, any 
pair of adjacent rods always has opposite polarity (Fig. 3). 
The potential inside the trap near the central axis (i.e., 
x, y <C R) is given by 


V{x } y) 


[V 0 (x 2 - y 2 ) cos (Qt)\ 
2 R 2 


(1) 


where x and y denote the particle’s position in the radial 
direction, R represents the distance between the central 
axis and the surface of the rods, and ft is the driving fre- 
quency. The resultant time-varying electric fields produce 
a trapping force that increases linearly in the radial direc- 
tion. A single charged particle within the trapping region 
will be subjected to a time-a veraged f orce towards the cen- 
tral axis of the trap (r = sjx 1 + y 2 = 0) where the field 
intensity is a minimum. A nonzero time-averaged restoring 
force for a single particle is derivable from a pseudopoten- 
tial, the effective potential existing in an inhomogeneous 
electric rf field [4] 
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In the case of an infinitely long linear trap, Eq. (2) modifies 
to 
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where q is the charge on the particle and x and y are the 
particle’s position in the radial direction averaged over one 
period T = (2tt/Q). The force generated from this pseu- 
dopotential produces particle dynamics that superimpose 
a slower oscillating motion characterized by frequency w, 
upon a faster micromotion due to the driving frequency Q. 
Stable confinement of the particles occurs when w < fi. 
The frequency w for radial particle oscillation can be de- 
rived from Eq. (3) as 



Although the particle is in dynamic equilibrium along the 
radial direction, it is still free to traverse the longitudinal 
( 2 ) axis. If there is no constraining force along this axis 
and the trap contains more than one particle, the coulomb 
repulsion between the particles tends to propel them out 
the ends of the trap. To counterbalance these forces, dc- 
biased electrodes are incorporated at either end of the trap 
in both ion and macroparticle systems. 


lenses or mirrors. The pseudopotential arising from the 
rf fields and the potential generated from the endcaps are 
solved separately. Superimposing these two solutions pro- 
duces the net potential everywhere inside the cylindrical 
trap. 

The pseudopotential is determined by employing the 
four-sectored cylindrical model in Fig. 2 with grounded 
endcaps and a surface voltage V 0 applied such that V o (0) = 
— V o {0 + tt/ 2). This voltage arrangement represents the 
applied voltage on the four-rod quadrupole at a given in- 
stant. Starting with Laplace’s equation in cylindrical coor- 
dinates, a general solution of the product form is assumed: 

<Hr,0,z) = R(r)Q(0)Z(z) (5) 

It can be shown that for this type of geometry [5] 
z) — [A co s(mO) -f B sin(r7?0)][C/ m (kr) 


III. Analytical Model for the Potential Inside 
a Linear Trap 

The complete electric field inside the linear ion trap 
consists of a static field arising from applied voltage on 
the endcaps and the time-varying fields resulting from a 
voltage arrangement on the trap rods. Unfortunately, the 
analytic determination of the net electric field for a finite- 
length four- rod linear trap is formidable and, a s yet, only 
an approximation of the field exists by considering the trap 
to be infinitely long. This approximation neglects the ef- 
fect trap ends have on the generated rf fields that influence 
the dynamics of trapped ions. Rather than neglecting the 
end effects, in the model described here the fields generated 
by the four-rod geometry are approximated by employing a 
hollow cylinder partitioned into four sections with endcaps 
at either end (Fig. 2). While providing a rough approxima- 
tion to the four-rod trap, the model primarily represents a 
new linear-trap geometry that can be accurately modeled 
to determine the exact confining field, including any signif- 
icant end effects. An exact analytic determination of the 
fields is highly useful since the shape of the trapped ion 
cloud depends strongly on the fields surrounding it. Be- 
cause the ions contribute significantly to the electric field 
within the trap, only certain density profiles are accept- 
able. With a complete picture of the fields generated by 
the cylindrical shell and endcaps, Poisson’s equation can 
be solved numerically to investigate the shape of the ion 
plasma. By virtue of its design, the cylindrical geometry 
also conveniently provides shielding against charge buildup 
that might occur in regions external to the trap such as 


+ DK m (lcr)][Esin(kz) + Fcos(kz)] (6) 

where I m and I\ m are the modified Bessel functions. The 
remaining coefficients are constants to be determined. 
With the requirements that F, D = 0 and 


the equation for $ takes on the modified general form 

oo oo 

<P(r,0,z)= ^ F [(4n cosjinO) 

m — 0 n nO 

+ D mn sin(m0)) ^/ m (—7-)) 

M™))] <•> 

The asymmetric voltage is applied to the four cylinder 
sections on the cylindrical surface (r — R) such that: 

f V 0 , -7T<0<^) 


-v 0 . 

=£- < 0 < 0 

Vo, 

0 < 0 < I 

-V 0 , 

% <0 < 7T 


The potential on the surface is thus given as 
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C(r,9,z) = 
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with the coefficients 
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Substituting these constants back into Eq. (8), and mul- 
tiplying $ by cos(Qf) to take into account that it varies 
with a frequency the time-varying potential due to the 
surface voltage is determined: 
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where m = 2, 6, 10, . . . and n = 1,3,5,.... Providing the 
secular frequency is much less than the frequency of mi- 
cromotion (cj Q), the macromotion of the particle can 
be described by the time-independent pseudopotential dis- 
cussed in Section II 


<i31 

and incorporating Eq. (12), 

ip — 4 ~q 2 (A 3 (r,0, z) + B 2 (r,0,z) + C 2 (r,0,z)) (14) 


where 
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(15) 


The inhomogeneous electric field arising from this effective 
potential produces a net force on a single particle that, 
averaged over one period of the micromotion (T = Q/ 27 t), 
is directed towards r — 0 in the radial direction and out the 
ends of the trap in the z direction. The z dependence of the 
pseudopotential is demonstrated by Figs. 4(a) and 4(b), 
which represent examples of at two different positions 
of 0. In both of the examples, = R/S with V ac = 200 V, 
Q/2ir = 60 Hz, R — 0.008 m, and L = 0.06 m. A 5 /im- 
diameter alumina particle with an approximate positive 
charge-to-mass ratio of 0.0053 C/kg is used to evaluate 
Eq. (14). 

The magnitude of t}>( 0 ) for fixed values of r and z is 
sinusoidal and oscillates from a maximum in regions near 
the boundaries of each of the four sections of the cylinder 
(0 = 0, 7r/2 , 7 t,2t t) to a minimum along the centers of the 
cylinder sections (0 — tt/4, 37t/4, 5tt /4, T 7t/ 4) (Fig. 5). It is 
interesting to note that the pseudopotential in the cylin- 
drical model indicates an electric field in 0 that is actually 
45 deg out of phase with the electric field generated in the 
linear four-rod trap. In the latter case, the pseudopoten- 
tial is a maximum in the vicinity of the rods (located at 
0 = 7t/4, 37t/4, 5tt/4 } 77t/ 4) and a minimum between them. 
This discrepancy arises because the spatial separation of 
the rods results in weaker fields between any two rods in 
contrast to the fields between the cylinder sections whose 
boundaries meet and give rise to large electric fields. The 
fact that the pseudopotential is a minimum at the section 
centers in the cylindrical linear trap suggests that small 
viewing holes may be added in this vicinity without sig- 
nificant disturbance of the fields within. 


In order to obtain the total potential inside the cylin- 
der, the potential arising from the endcaps must be deter- 
mined and superimposed on the pseudopotential. This is 
accomplished by solving Laplace’s equation for a cylinder 
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with a grounded surface and endcaps held at a constant 
voltage V e . The resulting potential inside the cylinder for 
this configuration is 

Jo(^) ( cosb ( Pm( -JT^) ) 

PmM^n) [ COsh(%£) j 


oo 

= 2V e J2 

m = 1 


(16) 


potential due to the dc ofTset is solved in a similar manner 
as the problem in Section II, and the results are 


4>(r,0,z) 
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n odd 
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4 Mtt) 


sin(m0) 


(18) 


where Jq and J\ are the zeroth and first Bessel func- 
tions, respectively, and P m is the mth root of Jo- Adding 
Eqs. (14) and (16) determines the net potential inside a 
finite cylinder of length L and radius R 

Vnet(r,0,z) = 'iP + $ e (17) 

The net potential can be interpreted geometrically by 
observing that it is composed of the pseudopotential whose 
curve, as demonstrated in Figs. 4(a) and 4(b), rises for 
fixed values of r and 0 from a minimum at z = 0 and 
z = L to a maximum in the center of the trap, and an 
“cndcap” potential that varies from a constant \\ on the 
ends to a minimum in the center of the trap (Fig. 6). By 
controlling the amount of voltage applied to the endcaps, 
xj/ net can be adjusted. The model suggests that there is an 
optimum endcap voltage that maximizes the continuity of 
the net potential in the z direction. Figure 7 demonstrates 
the z dependence of ^ n et when $ is greater than ip. 


IV. Analytical Model for a DC Field 
Perpendicular to the Cylindrical Axis 

In a horizontal configuration, the introduction of a 
dc voltage on the top two rods of the linear quadruple 
trap generally serves to counter the forces of gravity on 
macroparticles. Because the rods are of finite length, how- 
ever, the electric field is not directed exclusively in the ra- 
dial direction as in the case of infinite rod length. Instead, 
a distortion of the field lines occurs at the ends of the rods, 
giving rise to a field component in the longitudinal direc- 
tion. The effects on particle dynamics of an asymmetiic 
dc bias may be understood by again considering the four- 
rod quadrupole to be a hollow cylindrical trap split into 
two sections; the z axis is coincident with the axis of the 
linear quadrupole trap. Figure 8 depicts the cylindrical 
model with radius R and length L. To model the effects 
of the ofTset in the trap, a dc voltage is applied on the top 
half of the cylinder and the bottom half and endcaps are 
grounded. This potential can be added to developed 
in the previous section to provide a full solution of the po- 
tential inside a linear trap when a dc bias is added. The 


Figure 9 represents a computer evaluation for this func- 
tion of the first 22 values of m and n, respectively. Near 
the central axis of the trap (r = 0), the function varies 
smoothly, falling off to zero at z = 0 and z = L. The 
resultant electric field at r = 0 is 


dz 


-2V 0 

L 


oo 



n odd 



The surface voltage V 0 is inherently positive in Eq. (19). 
By referring to Figs. 10(a) and (b), it is evident that if a 
single negatively charged particle is subjected to this field, 
the force on that particle will be directed inward towards 
the center of the trap (z = L/ 2), growing abruptly from 
zero to a maximum as it approaches either end. The poten- 
tial inside the trap could have been determined similarly 
by solving Laplace's equation for a voltage V 0 applied to 
the bottom half of the cylinder instead of the top. In this 
case, if the particle had been taken to be positive, the net 
effect would still support the particle against gravity, but 
would accelerate the particle out one of the ends of the 
cylinder. 

The analogy to the linear cylindrical trap is apparent. 
In order to counterbalance the gravitational effects on the 
particle, the dc voltage applied to the top rods must be of 
the opposite sign to that of the particle. But as demon- 
strated above, this condition inherently gives rise to an 
inhomogeneous force on the particle that is directed longi- 
tudinally inward, towards the center of the trap. Likewise, 
a gravitational offset can also be accomplished by apply- 
ing a voltage to the bottom two trap rods that is the same 
polarity as the particle. This, however, additionally pro- 
duces an electric field that acts to force the particle out 
the ends of the trap. 

Farther out along the radial direction, the potential in- 
creases in magnitude and the function assumes a more 
rectangular shape. In the limit at r = R , V 0 spans the 
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rectangular shape. In the limit at r = R } V 0 spans the 
length of the cylinder and the function discontinuously 
jumps to 0 at either end. The electric field evaluated on 
the surface is 



where 6 is again taken to be a constant at tt/ 2. The model 
indicates that the electric field in the z direction towards 
the center of the trap depends on the radial position of the 
particle. This dependence is demonstrated by comparing 
Figs. 10(a) and (b). 

The cylindrical model approximates the field accurately 
and yields insight into the influences on the particle in the 
longitudinal direction from the gravitational ofFset. For 
the case of many trapped particles, the model suggests 
that in some instances a bias on the end electrodes is not 
really necessary to counterbalance coulomb repulsion be- 
tween particles, if a dc offset on the top rods is added. In 
this case, the particles will distribute themselves in such 
a manner as to create an equilibrium situation, balancing 
the coulomb forces and the force arising from the applied 
voltage. 


V. Conclusion 

This article has presented a model for the rf fields in a 
linear trap by considering a four-sector cylindrical geome- 
try. This geometry approximates a finite linear trap with 


rod electrodes. It also represents a new configuration for 
an ion trap that is amenable to analytic determination of 
the confining fields. Using the model for the finite-length 
linear cylindrical trap, it was shown that confining fields 
have a longitudinal component as well as a radial compo- 
nent. The z component of the rf electric field produces 
an extra term in the pseudopotential not present in the 
ideal, infinitely long trap. The details of the end effects 
are important in ion confinement. The analytic model in- 
dicates that the pseudopotential arising from the rf fields 
gives rise to a time-averaged force towards the center of 
the trap in the radial direction and away from the center 
in the 2 direction. This latter effect, combined with the 
ions’ own coulomb interaction, will accelerate the parti- 
cles out of the trap if biased endcaps are not present. The 
density of the ion cloud can be determined numerically by 
solving Poisson’s equation for the fields generated by the 
trap electrodes encompassing the ion plasma. Future re- 
search is planned for the evaluation of the plasma shape 
in a linear cylindrical trap. 

When a dc bias potential is applied to the trap elec- 
trodes (for example, in the case of a macroparticle trap 
where the gravitational force is counterbalanced with a dc 
offset), an electric field is generated in both the radial and 
longitudinal directions. An analytic model has been de- 
veloped that utilizes a cylindrical geometry to determine 
these fields. According to the model, a particle of the op- 
posite sign to an applied dc voltage on the upper electrodes 
will feel an inwardly directed force along the 2 axis in addi- 
tion to the transverse force. On the other hand, a particle 
with the same polarity as an applied dc voltage on the 
lower electrodes will be forced upwards in the transverse 
direction and towards the trap ends in the longitudinal 
direction. 
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Fig. 4. The pseudopotential ^(z) at two positions of : (a) z de- 
pendence of xj)(rfi % z) evaluated at r— fl/8 and 0 — 0 for a single 
alumina particle with e/m = 0.0053 C/kg and trapping parameters 
V ac = 200 V and Ql2ir = 60 Hz; and (b) \p(r,0,z) evaluated at r = 

Rl 8 and 0 = 7T/4. 



Fig. 5. Angular dependence of il'(r t O t z) evaluated at RIA and 
z — LI 2 for a single alumina particle with e/m = 0.0053 C/kg and 
trapping parameters V ac = 200 V and fi/2T — 60 Hz. 



Fig. 6. The longitudinal (z) dependence of $ e (r,z) evaluated at 
r= 0 with end electrode voltage V n = 20 V. 
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Fig. 10. The electric field E(z) within the cylindrical model: (a) at 
a fixed r= 0, 0 = tt/ 2, and V D = 20 V; and (b) at a fixed r= 3/7/4, 
0 = 7T/2, and V Q = 20 V. 





